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We consider a new quantum gate mechanism based on electron spins in coupled semiconductor 
quantum dots. Such gates provide a general source of spin entanglement and can be used for 
quantum computers. We determine the exchange coupling J in the effective Heisenberg model as 
a function of magnetic (B) and electric fields, and of the inter-dot distance a within the Heftier- 
London approximation of molecular physics. This result is refined by using sp-hybridization, and by 
the Hund-Mulliken molecular-orbit approach which leads to an extended Hubbard description for 
the two-dot system that shows a remarkable dependence on B and a due to the long-range Coulomb 
interaction. We find that the exchange J changes sign at a finite field (leading to a pronounced jump 
in the magnetization) and then decays exponentially. The magnetization and the spin susceptibilities 
of the coupled dots are calculated. We show that the dephasing due to nuclear spins in GaAs can 
be strongly suppressed by dynamical nuclear spin polarization and/or by magnetic fields. 
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I. INTRODUCTION 

Semiconductor quantum dots, sometimes referred to as 
artificial atoms, are small devices in which charge carri- 
ers are confined in all three dimensions jjj. The confine- 
ment is usually achieved by electrical gating and/or etch- 
ing techniques applied e.g. to a two-dimensional electron 
gas (2DEG). Since the dimensions of quantum dots are 
on the order of the Fermi wavelength, their electronic 
spectrum consists of discrete energy levels which have 
been studied in great detail in conductance and spec- 
troscopy measurements jljlflj] • In GaAs heterostructures 
the number of electrons in the dots can be changed one 
by one starting from zero ||. Typical laboratory mag- 
netic fields (B w IT) correspond to magnetic lengths 
on the order of Ib ~ 10 nm, being much larger than the 
Bohr radius of real atoms but of the same size as artifi- 
cial atoms. As a consequence, the dot spectrum depends 
strongly on the applied magnetic field |3J-Q]- In coupled 
quantum dots which can be considered to some extent 
as artificial molecules, Coulomb blockade effects || and 
magnetization have been observed as well as the for- 
mation of a delocalized "molecular state" || . 

Motivated by the rapid down-scaling of integrated cir- 
cuits, there has been continued interest in classical logic 
devices made of electrostatically coupled quantum dots 
H . More recently, the discovery of new principles of com- 
putation based on quantum mechanics |lC| ] has led to the 
idea of using coupled quantum dots for quantum compu- 
tation | |nj ; many other proposed implementations have 
been explored, involving NMR Ji~2|^|T4 1 , trapped ions [jl5| , 
cavity QED p6| , and Josephson junctions [jl7|. Solid- 
state devices open up the possibility of fabricating large 
integrated networks which would be required for realistic 
applications of quantum computers. A basic feature of 
the quantum-dot scenario [fllf is to consider the electron 



spin S as the qubit (the qubit being the basic unit of 
information in the quantum computer). This stands in 
contrast to alternative proposals also based on quantum 
dots -^lj , in which it is the charge (orbital) degrees of 
freedom out of which a qubit is formed and represented 
in terms of a pseudospin-1/2. However, there are two im- 
mediate advantages of real spin over pseudospin: First, 
the qubit represented by a real spin- 1/2 is always a well 
defined qubit; the two-dimensional Hilbcrt space is the 
entire space available, thus there are no extra dimensions 
into which the qubit state could "leak" J2^] . Second, dur- 
ing a quantum computation phase coherence of the qubits 
must be preserved. It is thus an essential advantage of 
real spins that their dephasing times in GaAs can be on 
the order of microseconds |2j|, whereas for charge de- 
grees of freedom dephasing times are typically much less, 
on the order of nanoseconds p^j25[ |. 

In addition to a well defined qubit, we also need a con- 
trollable "source of entanglement" , i.e. a mechanism by 
which two specified qubits at a time can be entangled 
p6[ so as to produce the fundamental quantum XOR (or 
controllcd-NOT) gate operation, represented by a uni- 
tary operator Uxor This can be achieved by tem- 
porarily coupling two spins fllf| . As we will show in de- 
tail below, due to the Coulomb interaction and the Pauli 
exclusion principle the ground state of two coupled elec- 
trons is a spin singlet, i.e. a highly entangled spin state. 
This physical picture translates into an exchange cou- 
pling J(t) between the two spins Si and S2 described by 
a Heisenberg Hamiltonian 



H s {t) = J(t) Si ■ S 2 . 



(1) 



If the exchange coupling is pulsed such that / dtJ(t)/h = 
Jot s /Ti = it (mod 2tt), the associated unitary time evo- 
lution U(t) = Texp(i /' H s (t)(It /%) corresponds to the 
"swap" operator U sw which simply exchanges the quan- 
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turn states of qubit 1 and 2 111). Furthermore, the quan- 
tum XOR can be obtained [[Ll| by applying the sequence 

exp(i(7r/2)5f ) exp(-i( 7 r/2)S*|)[/ s 1 i 2 exp(i7r Sl)Ull 2 = 

1/2 

Uxor, he. a combination of "square-root of swap" U S w 
and single-qubit rotations expiinSf), etc. Since Uxor 
(combined with single-qubit rotations) is proven to be a 
universal quantum gate |l^,|26|], it can therefore be used 
to assemble any quantum algorithm. Thus, the study of 
a quantum XOR gate is essentially reduced to the study 
of the exchange mechanism and how the exchange cou- 
pling J{t) can be controlled experimentally. We wish to 
emphasize that the switchable coupling mechanism de- 
scribed in the following need not be confined to quan- 
tum dots: the same principle can be applied to other 
systems, e.g. coupled atoms in a Bravais lattice, over- 
lapping shallow donors in semiconductors such as P in Si 
|28| , and so on. The main reason to concentrate here on 
quantum dots is that these systems are at the center of 
many ongoing experimental investigations in mesoscopic 
physics, and thus there seems to be reasonable hope that 
these systems can be made into quantum gates function- 
ing along the lines proposed here. 

In view of this motivation we study in the following the 
spin dynamics of two laterally coupled quantum dots con- 
taining a single electron each. We show that the exchange 
coupling J(B, E, a) can be controlled by a magnetic field 
B (leading to wave function compression) , or by an elec- 
tric field E (leading to level detuning), or by varying 
the barrier height or equivalently the inter-dot distance 
2a (leading to a suppression of tunneling between the 
dots). The dependence on these parameters is of direct 
practical interest, since it opens the door to tailoring the 
exchange Jit) for the specific purpose of creating quan- 
tum gates. We further calculate the static and dynamical 
magnetization responses in the presence of perpendicular 
and parallel magnetic fields, and show that they give ex- 
perimentally accessible information about the exchange 
J. Our analysis is based on an adaptation of Heitlcr- 
London and Hund-Mulliken variational techniques |2£| 
to parabolically confined coupled quantum dots. In par- 
ticular, we present an extension of the Hubbard approx- 
imation induced by the long-range Coulomb interaction. 
We find a striking dependence of the Hubbard parame- 
ters on the magnetic field and inter-dot distance which is 
of relevance also for atomic-scale Hubbard physics in the 
presence of long-range Coulomb interactions. Finally, we 
discuss the effects of dephasing induced by nuclear spins 
in GaAs and show that dephasing can be strongly re- 
duced by dynamically polarizing the nuclear spins and/or 
by magnetic fields. 

The paper is organized as follows. In Sec. II we intro- 
duce the model for the quantum gate in terms of coupled 
dots. In Sec. Ill we calculate the exchange coupling first 
in the Heitler-London and then in the Hund-Mulliken 
approach. There we also discuss the Hubbard limit and 
the new features arising from the long range nature of 
the Coulomb interactions. In Sec. IV we consider the 



effects of imperfections leading to dephasing and gate er- 
rors; in particular, we consider dephasing resulting from 
nuclear spins in GaAs. Implications for experiments on 
magnetization and spin susceptibilities are presented in 
Sec. V, and Sec. VI contains some concluding remarks 
on networks of gates with some suggestions for single- 
qubit gates operated by local magnetic fields. Finally, 
we mention that a preliminary account of some of the 
results presented here has been given in Ref. p0[. 
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FIG. 1. Two coupled quantum dots with one valence elec- 
tron per dot. Each electron is confined to the xy plane. The 
spins of the electrons in dots 1 and 2 are denoted by Si and 
S2. The magnetic field B is perpendicular to the plane, i.e. 
along the z axis, and the electric field E is in-plane and along 
the x axis. The quartic potential is given in Eq. and is 
used to model the coupling of two harmonic wells centered at 
(±0, 0, 0). The exchange coupling J between the spins is a 
function of B, E, and the inter-dot distance 2a. 

II. MODEL FOR THE QUANTUM GATE 

We consider a system of two laterally coupled quantum 
dots containing one (conduction band) electron each, see 
Fig. [I]. It is essential that the electrons are allowed to 
tunnel between the dots, and that the total wave func- 
tion of the coupled system must be antisymmetric. It is 
this fact which introduces correlations between the spins 
via the charge (orbital) degrees of freedom. For defi- 
niteness we shall use in the following the parameter val- 
ues recently determined for single GaAs heterostructure 
quantum dots || that are formed in a 2DEG; this choice 
is not crucial for the following analysis but it allows us 
to illustrate our analytical results with realistic numbers. 
The Hamiltonian for the coupled system is then given by 

H = ^ + C + H z = H mh + H z , 

i=l,2 

hi = ^- (pi- -A(ri)) +ex i E + V(r i ), (2) 
2m \ c J 

e 2 

C=- - 

k |ri - r 2 | 

The single-particle Hamiltonian hi describes the electron 
dynamics confined to the xy-pl&ne. The electrons have 
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an effective mass m (m = 0.067 m e in GaAs) and carry 
a spin-1/2 S^. The dielectric constant in GaAs is k = 
13.1. We allow for a magnetic field B = (0, 0, B) applied 
along the z-axis and which couples to the electron charge 
via the vector potential A(r) = ^(—y,x,0). We also 
allow for an electric field E applied in-plane along the x- 
direction, i.e. along the line connecting the centers of the 
dots. The coupling of the dots (which includes tunneling) 
is modeled by a quartic potential, 



V(x,y) = 



y 



(3) 



which separates (for x around ±a) into two harmonic 
wells of frequency ujq, one for each dot, in the limit of 
large inter-dot distance, i.e. for 2a 3> 2<ib, where a is 
half the distanc e between the centers of the dots, and 
<2b = y/h/muJo is the effective Bohr radius of a single 
isolated harmonic well. This choice for the potential is 
motivated by the experimental fact that the spectrum 
of single dots in GaAs is well described by a parabolic 
confinement potential, e.g. with hto Q — 3meV f|. We 
note that increasing (decreasing) the inter-dot distance 
is physically equivalent to raising (lowering) the inter-dot 
barrier, which can be achieved experimentally by e.g. ap- 
plying a gate voltage between the dots Thus, the ef- 
fect of such gate voltages is described in our model simply 
by a change of the inter-dot distance 2a. We also note 
that it is only for simplicity that we choose the two dots 
to be exactly identical, no qualitative changes will occur 
in the following analysis if the dots are only approxi- 
mately equal and approximately of parabolic shape. 

The (bare) Coulomb interaction between the two elec- 
trons is described by C. The screening length A in almost 
depleted regions like few-electron quantum dots can be 
expected to be much larger than the bulk 2DEG screen- 
ing length (which is about 40 nm in GaAs). Therefore, 
A is large compared to the size of the coupled system, 
A ^> 2a w 40 nm for small dots, and we will consider 
the limit of unscreened Coulomb interaction (A/a S> 1) 
throughout this work. 

The magnetic field B also couples to the electron spins 
via the Zeeman term Hz = g/ib J2i B-t ■ Si, where g is 
the effective g-factor (g w —0.44 for GaAs), and /xb the 
Bohr magneton. The ratio between the Zeeman split- 
ting and the relevant orbital energies is small for all B- 
values of interest here; indeed, gn^B/TwJo < 0.03, for 
B <C Bq — (hoJo/ IMB){m/m e ) ~ 3.5T, and g[i^B /Hlul ^$ 
0.03, for B 3> Bq, where ojl = eB/2mc is the Lar- 
mor frequency, and where we used htoo = 3meV. Thus, 
we can safely ignore the Zeeman splitting when we dis- 
cuss the orbital degrees of freedom and include it later 
into the effective spin Hamiltonian. Also, in the few- 
electron system we are dealing with, spin-orbit effects can 
be completely neglected since H so /Uujo rs 10~ 7 , where 
H so = (wg/2mc 2 )L-S is the spin-orbit coupling of an elec- 
tron in a parabolic confinement potential pCf , This has 
the important implication that dephasing effects induced 



e.g. by potential or charge fluctuations in the surround- 
ings of the isolated dots can couple only to the charge of 
the electron so that they have very small influence on the 
phase coherence of the isolated spin itself (for dephasing 
induced by coupling the dots see Sec. |fv| ). It is for this 
reason that it is preferable to consider dots containing 
electrons instead of holes, since holes will typically have 
a sizable spin-orbit interaction [JjJ . 

Finally, we assume a low-temperature description 
where kT -C Tvjjq, so that we can restrict ourselves to 
the two lowest orbital eigenstates of H OI \>, one of which is 
symmetric (spin singlet) and the other one antisymmetric 
(spin triplet). In this reduced (four-dimensional) Hilbcrt 
space, -fforb can be replaced by the effective Heisenberg 
spin Hamiltonian Eq. (Q), H s = JSi ■ S2, where the ex- 
change energy J = e t — e s is the difference between the 
triplet and singlet energy which we wish to calculate. The 
above model cannot be solved in an analytically closed 
form. However, the analogy between atoms and quantum 
dots (artificial atoms) provides us with a powerful set of 
variational methods from molecular physics for finding e t 
and e s . Note that the typical energy scale hajQ ss meV in 
our quantum dot is about a thousand times smaller than 
the energies (Ry « eV) in a hydrogen atom, whereas the 
quantum dot is larger by about the same factor. This is 
important because their size makes quantum dots much 
more susceptible to magnetic fields than atoms. In anal- 
ogy to atomic physics, we call the size of the electron 
orbitals in a quantum dot the Bohr radius, although it is 
determined by the confining potential rather than by the 
Coulomb attraction to a positively charged nucleus. For 
harmonic confinement ob = \Jh/mu)Q is about 20 nm for 
hu>o = 3 meV. 



III. EXCHANGE ENERGY 

A. Heitler-London approach 

We consider first the Heitler-London approximation, 
and then refine this approach by including hybridiza- 
tion as well as double occupancy in a Hund-Mulliken 
approach, which will finally lead us to an extension of 
the Hubbard description. We will see, however, that 
the qualitative features of J as a function of the con- 
trol parameters are already captured by the simplest 
Heitler-London approximation for the artificial hydrogen 
molecule described by Eq. ^. In this approximation, one 
starts from single-dot ground-state orbital wavefunctions 
(p(r) and combines them into the (anti-) symmetric two- 
particle orbital state vector 



112) ±121) 
v/2(l±S 2 ) : 



(4) 



the positive (negative) sign corresponding to the spin sin- 
glet (triplet) state, and S = j d 2 rtp* +a (r)ip^ a (r) = (2|1) 
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denoting the overlap of the right and left orbitals. A non- 
vanishing overlap implies t hat th e electrons tunnel be- 
tween the dots (see also Sec. [II B ). Here, tp_ a (r) = (rjl) 
and <^+ a (r) = (r|2) denote the one-particle orbitals cen- 
tered at r = (=pa, 0), and \ij) = \i)\j) are two-particle 
product states. The exchange energy is then obtained 
through J = e t — e s = (*_ |ff rb|*-) - (V + \H OTb \V + ). 
The single-dot orbitals for harmonic confinement in two 
dimensions in a perpendicular magnetic field are the 
Fock-Darwin states |3l|| , which are the usual harmonic 
oscillator states, magnetica lly compressed by a factor 
6 = uj/luq = yl + uj 2 /u> 2 , where wl = eB/2mc de- 
notes the Larmor frequency. The ground state (energy 
huj = bhu>o) centered at the origin is 



J/meV 




2 4 6 



8 B/T 



muJ - m uj(x 2 +y 2 )/2h 



itTi 



(•5) 



Shifting the single particle orbitals to (±a, 0) in the 
presence of a magnetic field we obtain if± a {x 1 y) — 
exp(±iya/2lg)ip(x =p a,y)- The phase factor involv- 
ing the magnetic length lg = ^hc/eB is due to the 
gauge transformation A± a = B(—y,x =F a, 0)/2 — > A = 
B(— y, x, 0)/2. The matrix elements of i? rb needed to 
calculate J are found by adding and subtracting the 
harmonic potential centered at x = — (+) a for elec- 
tron 1(2) in iJorb, which then takes the form H ol -h = 
h°- a (rx) + h° +a (r 2 ) + W + C, where h° ±a ( ri ) = ( Pi - 
eA(ri)/c) 2 /2m + mw 2 ((i, q= a) 2 + yf)/2 is the Fock- 
Darwin Hamiltonian shifted to (±a,0), and W(x,y) — 
V(x, y) - mw 2 ((n + a) 2 + (x 2 - a) 2 )/2. We obtain 



J 



2S 2 



l-S 4 



12|C + W|12) 



Re(12|C + W\21) 

S 2 



(G) 



where the overlap becomes S = exp(— mwa 2 /% — 
a 2 h/ 'AlginLj). Evaluation of the matrix elements of C 
and W yields (see also pCf) 



J = 



sinh(2d 2 (26- 1/6)) 



cVbi e- bd2 l (bd 2 ) 



^ b -V b h (d 2 (b-l/b)) \ +±(l + b<?) 



(7) 



where we introduce the dimensionless distance d = o/ae, 
and I is the zeroth order Bessel function. The first and 
second terms in Eq. (fjj) are due to the Coulomb inter- 
action C, where the exchange term enters with a minus 
sign. The parameter c = ^ it / 2(e 2 / na^) / Tiujq (~ 2.4, for 
Hluq = 3 meV) is the ratio between Coulomb and confin- 
ing energy. The last term comes from the confinement 
potential W. The result J{B) is plotted in Fig. || (dashed 
line). Note that typically |J/7iwo| < 0.2. Also, we see 
that J > for B = 0, which must be the case for a 
two-particle system that is time-reversal invariant |29[| . 



FIG. 2. Exchange energy J in units of meV plotted 
against the magnetic field B (in units of Tesla), as obtained 
from the s-wave Heitler-London approximation (dashed line), 
Eq. (^), and the result from the improved sp-hybridized 
Heitler-London approximation (triangles) which is obtained 
numerically as explained in the text. Note that the qual- 
itative behavior of the two curves is similar, i.e. they 
both have zeroes, the s-wave approximation at Bl, and the 
sp-hybridized approximation at Bt P , and also both curves 
vanish exponentially for large fields. Bo = (fkJo/HB)('m/m e ) 
denotes the crossover field to magnetically dominated confin- 
ing (B 3> -Bo). The curves are given for a confinement energy 
Tiujo = 3meV (implying for the Coulomb parameter c — 2.42), 
and inter-dot distance o = 0.7 Ob- 

The most remarkable feature of J(B), however, is the 
change of sign from positive to negative at B = Bl, 
which occurs over a wide range of parameters c and a. 
This singlet-triplet crossing occurs at about Bl — 1.3 T 
for huj = 3meV (c = 2.42) and d = 0.7. The tran- 
sition from antiferromagnetic (J > 0) to ferromagnetic 
(J < 0) spin-spin coupling with increasing magnetic field 
is caused by the long-range Coulomb interaction, in par- 
ticular by the negative exchange term, the second term 
in Eq. (g). As B » B (~ 3.5 T for huj = 3meV), 
the magnetic field compresses the orbits by a factor 
6 ~ B/Bo > 1 and thereby reduces the overlap of 
the wavefunctions, S 2 « exp(— 2d 2 (2b — 1/6)), exponen- 
tially strongly. Similarly, the overlap decays exponen- 
tially for large inter-dot distances, d ^> 1. Note how- 
ever, that this exponential suppression is partly com- 
pensated by the exponentially growing exchange term 
(12\C\21)/S 2 cx exp(2d 2 (6 - 1/6)). As a result, the ex- 
change coupling J decays exponentially as exp(— 2<i 2 &) 
for large 6 or d, as shown in Fig. ^3 for B — (6=1). 
Thus, the exchange coupling J can be tuned through zero 
and then suppressed to zero by a magnetic field in a very 
efficient way. We note that our Heitler-London approxi- 
mation breaks down explicitly (i.e. J becomes negative 
even when B = 0) for certain inter-dot distances when c 
exceeds 2.8. Finally, a similar singlet-triplet crossing as 
function of the magnetic field has been found in single 
dots with two electrons fB2|. 
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The exchange energy J also depends on the applied 
electric field E. The additional term e[x\ + x^jE in 
the potential merely shifts the one-particle orbitals by 
Ax = eE/majQ, raising the energy of both the singlet 
and triplet states. Since the singlet energy turns out to 
be less affected by this shift than the triplet, the exchange 
energy J increases with increasing E, 

J(B, E) = J(B, 0) + smh(2 / ( 2° _ 1/b)) ~ (^) , 

(8) 

the increase being proportional to itiujq (Ax) 2 . [We note 
that this increase of J(B,E) is qualitatively consistent 
with what one finds from a standard two-level approxi- 
mation of a ID double-well potential (with J(B, 0) be- 
ing the effective tunnel splitting) in the presence of a 
bias given by eEa.] The variational Ansatz leading 
to Eq. (||) is expected to remain accurate as long as 
J(B,E)-J(B,0) < J(B,0); for larger ^-fields the levels 
of the dots get completely detuned and the overlap of the 
wavefunctions (i.e. the coherent tunneling) between the 
dots is suppressed. Of course, a sufficiently large elec- 
tric field will eventually force both electrons on to the 
same dot, which is the case when eEa exceeds the on- 
site repulsion U(^> J(B,E — 0), see below). However, 
this situation, which would correspond to a quantum dot 
helium [ p^3[ , is not of interest in the present context. Con- 
versely, in case of dots of different size (or shape) where 
the energy levels need not be aligned a priori, an appro- 
priate electric field can be used to match the levels of the 
two dots, thus allowing coherent tunneling even in those 
systems. Recent conductance measurements || on cou- 
pled dots of different size (containing several electrons) 
with electrostatic tuning have revealed clear evidence for 
a delocalized molecular state. 

A shortcoming of the simple approximation described 
above is that solely ground-state single-particle orbitals 
were taken into account and mixing with excited one- 
particle states due to interaction is neglected. This 
approximation is self-consistent if J <C Ae, where Ae 
denotes the single-particle level separation between the 
ground state and the first excited state. We find |J/Ae| < 
0.25 at low fields B < 1.75 T, therefore J(B) is at least 
qualitatively correct in this regime. At higher fields 
| J/Ae| ss 1, indicating substantial mixing with higher or- 
bitals. An improved Hcitlcr-London variational Ansatz is 
obtained by introducing sp-hybridized single-dot orbitals 
(in analogy to molecular physics), i.e. <\> = cp s + aip px + 
iPfpy, where <p s — ip is the s-orbital introduced above, 

<Ppq = \J^muqeyi i p(—rriLur' 2 /2Ti)/'h, q — x,y, are the low- 
est two Fock-Darwin excited states (at zero field) with an- 
gular momentum \£\ = 1, and a and (3 are real variational 
parameters to be determined by minimization of the sin- 
glet and triplet energies e Sit (a,/3), which is done numer- 
ically. The ippq are chosen to be real, they are however 
not eigenstates of the single-particle Hamiltonian, which 



are ip px ± i(p py (with eigenenergy 2fuv ± Tilo-^). Note that 
while e Sit decrease only by ~ 1% due to hybridization, the 
relative variation of J = e t — e s can still be substantial. 
Nevertheless, the resulting exchange energy J sp (Fig. ||) 
is only quantitatively different from the pure s-wave re- 
sult J = J s , Eq. (|). At low fields, J s p < J s and the 
change of sign occurs already at about Bl p ~ 0.4 T < Bl. 
At high fields, J sp shows a much more pronounced decay 
as a function of B. 

Being a completely orbital effect, the exchange interac- 
tion between spins of course competes with the Zeeman 
coupling Hz of the spins to the magnetic field. In our 
case, however, the Zeeman energy Hz is small and ex- 
ceeds the exchange energy (polarizing the spins) only in a 
narrow window (about 0.1 T wide) around Bt p and again 
for high fields (B > 4 T). 




FIG. 3. The exchange coupling J obtained from 
Hund-Mulliken (full line), Eq. (|ll|), and from the extended 
Hubbard approximation (dashed line), Eq. (^). For compari- 
son, we also plot the usual Hubbard approximation where the 
long-range interaction term V is omitted, i.e. J = <kt%/Uy{ 
(dashed-dotted line). In (a), J is plotted as a function of the 
magnetic field B at fixed inter-dot distance (d = a/as = 0.7), 
and for c = 2.42, in (b) as a function of inter-dot distance 
d — a/a,B at zero field (B = 0), and again c = 2.42. For these 
parameter values, the s-wave Heitler- London J, Eq. (Q), and 
the Hund-Mulliken J (full line) are almost identical. 
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B. Hund-Mulliken approach and Hubbard Limit 

We turn now to the Hund-Mulliken method of molec- 
ular orbits ]2{| which extends the Heitler-London ap- 
proach by including also the two doubly occupied states, 
which both are spin singlets. This extends the orbital 
Hilbert space from two to four dimensions. First, the 
single particle states have to be orthonormalized, lead- 
ing to the states $ ±a = (ip± a - g^ a )l^\ - 2Sg + g 2 , 
where S again denotes the overlap of (p- a with <p +a and 
g = (1 — vl — S 2 )/S. Then, diagonalization of 



H or b = 2e 



I u 


X 


-V2t R 





X 


u 







-y/2ta 


-y/2ta 


v+ 





V o 








V- 



(9) 



in the space spanned by \&± a (ri,r 2 ) = ^± a {ri)^± a (r 2 ), 
<r ± ( ri ,r 2 ) = [* +0 (n)$_ a (r 2 ) ± $_ (n)$ +a (r 2 )]/>/2 
yields the eigenvalues e s ± = 2e + E/h/2 + V+ ± 
v/C/h/4 + 4^, e s0 = 2e + J7 H - 2X + (singlet), and 
e t = 2e + V_ (triplet), where the quantities 

e=(<S>±a\h Q ±a \<S>±a), 

tB=t-w= ($ ±a i4i$ T „) - (nM*i a )/V2, 

v = k_ - v+ = (* s _|c|# s _) - (* s + |c|* h ;), (io) 

x={*i a \c\*% a ), 

U U = U -V++ X 

= - (nicin> + (*lici^«>, 

all depend on the magnetic field B. The exchange energy 
is the gap between the lowest singlet and the triplet state 



J = e t 



U u , 1 
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(11) 



In the standard Hubbard approach for short-range 
Coulomb interactions (and without -B-field) p{| J re- 
duces to — U /2 + \/U 2 + 16i 2 /2, where f denotes the hop- 
ping matrix element, and U the on-site repulsion (cf. Eq. 
(|lO|)). Thus, £h and Uh are the extended hopping matrix 
element and the on-site repulsion, resp., renormalized by 
long-range Coulomb interactions. The remaining two sin- 
glet energies e s + and e s o are separated from e t and e s _ 
by a gap of order [/r and are therefore neglected for the 
study of low-energy properties. The evaluation of the 
matrix elements is straightforward but lengthy, and we 
give the results in Appendix |A|. Typically, the "Hubbard 
ratio" tu/Un is less than 1, e.g., if d = 0.7, Tiuj = 3meV, 
and B — 0, we obtain ta/Uu = 0.34, and it decreases 
with increasing B. Therefore, we are in an extended 
Hubbard limit, where J takes the form 



At 2 



(12) 



The first term has the form of the standard Hubbard ap- 
proximation |35[ (invoked previously \L% ) but with in 



and Uu being renormalized by long-range Coulomb in- 
teractions. The second term V is new and accounts for 
the difference in Coulomb energy between the singly oc- 
cupied singlet and triplet states . It is precisely this V 
that makes J negative for high magnetic fields, whereas 
t^/Un > for all values of B (see Fig. |a). Thus, the 
usual Hubbard approximation (i.e. without V) would 
not give reliable results, neither for the independence 
(Fig. |^a) nor for the dependence on the inter-dot distance 
a (Fig. ||b) j36|. Since only the singlet space has been en- 
larged, it is clear that we obtain a lower singlet energy 
e s than that from the s-wave Heitler-London calculation, 
but the same triplet energy et, and therefore J = et — e s 
exceeds the s-wave Heitler-London result, Eq. How- 
ever, the on-site Coulomb repulsion U oc c strongly sup- 
presses the doubly occupied states ^± a and already for 
the value of c = 2.4 (corresponding to Tiljq — 3meV) we 
obtain almost perfect agreement with the s-wave Heitler- 
London result (Fig. ||). For large fields, i.e. B ^> Bq, the 
suppression becomes even stronger (U oc VB) because 
the electron orbits become compressed with increasing 
B and two electrons on the same dot are confined to a 
smaller area leading to an increased Coulomb energy. 



IV. DEPHASING AND QUANTUM GATE 
ERRORS 

We allow now for imperfections and discuss first the 
dephasing resulting from coupling to the environment, 
and then address briefly the issue of errors during the 
quantum-gate operation. We have already pointed out 
that dephasing in the charge sector will have little ef- 
fect on the (uncoupled) spins due to the smallness of 
the spin-orbit interaction. Similarly, the dipolar interac- 
tion between the qubit spin and the surrounding spins is 
also minute, it can be estimated as (^Mb) 2 /^b ~ 10~ 9 
meV. Although both couplings arc extremely small they 
will eventually lead to dephasing for sufficiently long 
times. We have described such weak-coupling dephas- 
ing in terms of a reduced master equation elsewhere flil"| , 
and we refer the interested reader to this work. Since 
this type of dephasing is small it can be eliminated by 
error correction schemes f37f . 

Next, we consider the dephasing due to nuclear spins 
in GaAs semiconductors, where both Ga and As possess 
a nuclear spin / = 3/2. There is a sizable hyperfine 
coupling between the electron-spin (s = 1/2) and all the 
nuclear spins in the quantum dot which might easily lead 
to a flip of the electron spin and thus cause an error in the 
quantum computation. We shall now estimate this effect 
and show that it can be substantially reduced by spin 
polarization or by a field. We consider an electron spin 
S in contact with N nuclear spins I*- 1 -* in the presence of 
a magnetic field B \\ z. The corresponding Hamiltonian 
is given by H = AS ■ I + b z S z + b z I z = H + V , where 



G 



H = ASJ Z + b z S z + 6,7, , V = A(S+I- + S_7 + )/2. 

(13) 

Here, A is a hyperfine coupling, I = ^2i=i ^ ^ s the total 
nuclear spin, and b z = g^B z , b z = gN[iissB z (gpj and 
//jv denote the nuclear g factor and magneton). Con- 
sider the initial eigenstate \i) of Ho, which we will con- 
sider to be one basis vector for the qubit, where the 
electron spin is up (in the S z basis), and the nuclear 

spins are in a product state of T^-eigenstates with total 
I z = pNI (—1 < p < 1), i.e. in a state with polariza- 
tion p along the z-axis; here, p = ±1 means that the 
nuclear spins are fully polarized in positive (negative) 
z-direction, and p = means no polarization. Due to 
the hyperfine coupling the electron spin can flip (i.e. de- 
phase) with the entire system going into a final state 
which is again a product state but now with the electron- 
spin down, and, due to conservation of total spin, the 
^-component 7* of one and only one nuclear spin hav- 
ing increased by 2s = 1. All final states are degen- 
erate and again eigenstates of H with eigenenergy Ef. 
We will consider this process now within time-dependent 
perturbation theory and up to second order in V. The en- 
ergy difference between initial and final states amounts 
to Ei — Ef s» 2s[A(pIN + s) + b z ), where we use that 
b z 3> b z . For the reversed process with an electron- 
spin flip from down to up but with the same initial po- 
larization for the nuclear spins the energy difference is 
w — 2s[A(pIN — s) +6J. The total transition probability 
to leave the initial state \i) after time t has elapsed is 
then 



We interpret this total transition probability Pi (t) as the 
degree of decoherence caused by spin-flip processes over 
time t. Now, | (k\V \i) | 2 = A 2 [1(1 + 1) - li k) [li k) + l)]/4. 
Assuming some distribution of the nuclear spins we can 
replace this matrix element by its average value (denoted 



by brackets) where y ((T^) 2 ) describes then the vari- 
ance of the mean value (7j*^) = pi. E.g. a Poissonian dis- 
tribution gives |(fc|F|i)| 2 w A 2 [I(I+l)-pI(pI+l)}/4 : in 
which case the matrix element vanishes for full polariza- 
tion parallel to the electron-spin (i.e. p = 1), as required 
by conservation of total spin. Pi(t) is strongly suppressed 
for final states for which to = 2nh/\Ei — Ej\ -C t, which 
simply reflects conservation of energy. In particular, for a 
substantial nuclear polarization, i.e. p 2 N ^S> 1, Pi{t) os- 
cillates in time but with the vanishingly small amplitude 
1 jp 2 N (for B = 0) . We can estimate N to be on the order 
of the number of atoms per quantum dot, which is about 
10 5 . Such a situation with p 2 N ^> 1 can be established 
by dynamically spin-polarizing the nuclear spins (Over- 
hauser effect) e.g. via optical pumping |pq] or via spin- 
polarized currents at the edge of a 2DEG T39| . This gives 



rise to an effective nuclear field B n = ApNI /g/iB which is 
reported to be as large as B* = 4 T in GaAs (correspond- 
ing to p — 0.85) |3j| and which has a lifetime on the order 
of minutes [j38| . Alternatively, for unpolarized nuclei with 
p = but a field B in the Tesla range, the amplitude of 
Pi(t) vanishes as (AIN/ g^ B B) 2 /N m (B*/B) 2 /N < 1. 
For B or B n = 1 T the oscillation frequency 1/tg of Pi (t) 
is about 10 GHz. Thus, spin flip processes and hence de- 
phasing due to nuclear spins can be strongly suppressed, 
either by dynamically polarizing the nuclear spins and/or 
by applying a magnetic field B. The remaining dephasing 
effects (described again by a weak-coupling master equa- 
tion [[nj) should then be small enough to be eliminated 
by error correction. 

We now address the imperfections of the quantum gate 
operation. For this we note first that, for the purpose of 
quantum computing, the qubits must be coupled only 
for the short time of switching t s , while most of the time 
there is to be no coupling between the dots. We estimate 
now how small we can choose r s . For this we consider a 
scenario where J (initially zero) is adiabatically switched 
on and off again during the time r s , e.g. by an elec- 
trical gate by which we lower and then raise again the 
barrier V(t) between the dots (alternatively, we can vary 
B, a, or E). A typical frequency scale during switch- 
ing is given by the exchange energy (which results in the 
coherent tunneling between the dots) averaged over the 
time interval of switching, J = (l/r s ) J^ B dtJ(t). Adi- 
abaticity then requires that many coherent oscillations 
(characterized approximately by J) have to take place in 
the double-well system while the control parameter v = 
V, B, a, or E is being changed, i.e. 1/t s w \v/v\ <C J /h. 
If this criterion is met, we can use our equilibrium anal- 
ysis to calculate J(v) and then simply replace J(v) by 
J(v(t)) in case of a time-dependent control parameter 
v(t) |4l|]. Note that this is compatible with the require- 
ment needed for the XOR operation, Jr s /h — nn, n 
odd, if we choose n 3> 1. Our method of calculating 
J is self-consistent if J <§; Ae, where Ae denotes the 
single-particle level spacing. The combination of both 
inequalities yields l/r s <C J/h <C Ae/7i, i.e. no higher- 
lying levels can be excited during the switching. Finally, 
since typically J i=s 0.2 meV we see that r s should not 
be smaller than about 50 ps. Now, during the time r s 
spin and charge couple and thus dephasing in the charge 
sector described by can induce dephasing of spin via 
an uncontrolled fluctuation 8 J of the exchange coupling. 
However, this effect is again small, it can be estimated 
to be on the order of t s /t^ ~ 10~ 2 , since even for large 
dots is reported to be on the order of nanoseconds |24] . 
This seems to be a rather conservative estimate and one 
can expect the spin dephasing to be considerably smaller 
since not every charge dephasing event will affect the 
spin. Finally, weak dephasing of the effective spin Hamil- 
tonian during switching has been described elsewhere [Tl| ] 
in terms of a weak-coupling master equation which ac- 
counts explicitly for decoherence of the spins during the 
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switching process. Based on this analysis ]D| , the prob- 
ability for a gate error per gate operation (described by 
fC 2 in Eq. (13) of is estimated to be approximately 
Ts/t? ~ 10~ 2 or better (see above). 



J/meV 



M/|X E 




FIG. 4. The equilibrium magnetization M (box-shaped 
symbols) in units of Bohr magnetons /xb as a function of mag- 
netic field. M is obtained numerically from the sp-hybridized 
Heitler-London approximation. Note that the magnetization 
exhibits a jump at the field value Bl P for which the ex- 
change J sp (triangle symbols) changes sign. At the left and 
right hand side of the jump the negative slope of M(B) indi- 
cates orbital diamagnetism. The temperature for this plot is 
T = 0.2 K, while as before hu>o = 3meV and a = 0.7 as- 



V. EXPERIMENTAL IMPLICATIONS 

Coherent coupling between the states of neighboring 
dots is the keystone of our proposal for quantum gate op- 
eration, and experimental probes of this coupling will be 
very interesting to explore. The effect of the dot-dot cou- 
pling manifests itself in the level structure, which could 
be measured non-invasively with spectroscopic methods 
[pLpj. An alternative way is to measure the static mag- 
netization in response to a magnetic field B which is ap- 
plied along the z-axis. This equilibrium magnetization is 
given by M = gn B Tr(Sf + Si)e-( H ° +Hz V kT , where H s is 
given in Eq. (|]J), and Hz — g/iB Yli B< ■ Si is the Zeeman 
term. It is straightforward to evaluate M, and in Fig. |] 
we plot M as a function of B for a typical temperature 
T = 0.2 K. The exchange J sp (B) is also shown in Fig. |. 
Both J sp (B) and M are the results of the sp-hybridized 
Heitler-London approximation. We note that the equilib- 
rium magnetization M(B) is strongly dominated by the 
orbital response (via the exchange J) ; we find a diamag- 
netic response (negative slope of M) for B < Bt p which 
is followed by a pronounced jump in the magnetization at 
the field Bl p followed again by a diamagnetic response. 
Experimental observation of this jump would give ev- 
idence for the existence of the predicted singlet-triplet 
level-crossing at _B» P , and such measurements would al- 
low one to "map out" J around the point where it can be 
tuned to zero, e.g. by also varying the barrier between 



the dots. The magnetic moment produced by the orbital 
motion of the electrons in one pair of coupled quantum 
dots at the peak (B = Bl p ) is around 10/xb (see Fig. ||). 
This signal could be further amplified by using an ensem- 
ble of pairs of coupled quantum dots. 

A further way to get experimental information about 
the exchange coupling would be to measure the spin re- 
sponse to an ac magnetic field (in the linear-response 
regime), described by the dynamical spin susceptibilities 

X £?» = (t/ft) J™ dteMi"t)([SUt), S«(0)]), where 
m,n — 1,2, and p,q — x,y,z. Being interested in the 
spin response only, we assume this ac field to be applied in 
plane so that there is no orbital response (for a sufficiently 
weak field with no subband mixing). We see then that 
all the transverse spin susceptibilities Xnm' q vanish, and 
we are left with the longitudinal ones only, where Xmn = 
Xmn = Xmn = Xmn due to the rotational symmetry of H s . 
It is sufficient to consider the dissipative part, Xmni 10 ) = 
ImXmnM, for which we obtain Xii = X22 = ~Xi2 = 
-X 21 = -(n/4)f(J,B)[5(hu) + J) — 8{hw - J)], where 
f(J,B) = (e J / kT - 1)/(1 + e J ' kT + 2 cosh(g /ji B B/kT)). 
Also, due to conservation of total spin, the total response, 
Xij + X2j, as well as the response to a spatially uniform 
field, xn + Xi2j vanish. Thus, to observe the spin sus- 
ceptibilities calculated here one needs to apply the fields 
locally or to measure the spin of a dot separately; both 
cases could be realized e.g. by atomic or magnetic force 
microscopes (see also below, where we briefly discuss lo- 
cal fields produced by field gradients). 



VI. CONCLUDING REMARKS 

We end with a few comments on a network of cou- 
pled quantum dots in the presence of fields (see also Ref. 
[JTTJ] ) . In a set-up with only one quantum gate (i.e. two 
quantum dots) the gate operation can be performed us- 
ing uniform magnetic fields (besides electric gates) , while 
in a quantum computer with many gates, which have to 
be controlled individually, local magnetic fields are in- 
dispensable, especially for the single-qubit gates pr| . [42[ . 
However, we emphasize that it is not necessary that ev- 
ery single quantum dot in a network is directly address- 
able with a local magnetic field. Indeed, using "swap" 
operations U sw , any qubit-state can be transported to 
a region where the single-qubit gate operation is per- 
formed, and then back to its original location, without 
disturbing this or other qubits. In one possible mode 
of operation a constant field defined by J(-B*) = 0, 
is applied, while smaller time-dependent local fields then 
control the gate operations. We can envision local fields 
being achieved by a large number of techniques: with 
neighboring magnetic dots p| , closure domains, a grid 
of current-carrying wires below the dots, tips of magnetic 
or atomic force microscopes, or by bringing the qubit 
into contact (by shifting the dot via electrical gating) 
with a region containing magnetic moments or nuclear 
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spins with different hyperfine coupling (e.g. AlGaAs in- 
stead of GaAs)-and others. A related possibility would 
be to use magnetic field gradients. Single-qubit switch- 
ing times of the order of t s w 20 ps require a field of 

1 T, and for an inter-dot distance 2a w 30 nm, we would 
need gradients of about lT/30nm, which could be pro- 
duced with commercial disk reading/writing heads. [The 
operation of several XOR gates via magnetic fields also 
requires gradients of similar magnitude.] Alternatively, 
one could use an ac magnetic field _B ac and apply electron 
spin resonance (ESR) techniques to rotate spins with a 
single-qubit switching time (at resonance) t s « TtTi/ B ac . 
To address the dots of an array individually with ESR, 
a magnetic field gradient is needed which can be esti- 
mated as follows. Assuming a relative ESR linewidth of 
1% and again 2a = 30 nm we find about B ac ■ 10 4 cm _1 . 
Field gradients in excitation sequences for NMR up to 

2 ■ 10 4 G/cm have been generated which allows for 
_B ac « 1G. The resulting switching times, however, are 
rather long, on the order of 100 ns, and larger field gra- 
dients would be desirable. Finally, such ESR techniques 
could be employed to obtain information about the effec- 
tive exchange values J: the exchange coupling between 
the spins leads to a shift in the spin resonance frequency 
which we found to be of the order of J/h by numerical 
analysis Q. 

To conclude, we have calculated the exchange energy 
J(B, E, a) between spins of coupled quantum dots (con- 
taining one electron each) as a function of magnetic and 
electric fields and inter-dot distance using the Heitler- 
London, hybridized Heitler-London, and Hund-Mulliken 
variational approach. We have shown that J(B, E, a) 
changes sign (reflecting a singlet-triplet crossing) with in- 
creasing B field before it vanishes exponentially. Besides 
being of fundamental interest, this dependence opens up 
the possibility to use coupled quantum dots as quantum 
gate devices which can be operated by magnetic fields 
and/or electric gates (between the dots) to produce en- 
tanglement of qubits. 
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= eB/2mc. The single-particle matrix elements are 
given by 



3 1 



3 S 2 



1 



32 b 2 d 2 8 1 - S 2 \ b 



3 5/1 , 2 
t=-z — I T+d 2 



b, 



1 - S 2 \ b 



(Al) 
(A2) 



where we used S = exp(— c? 2 (26 — 1/6)). The (two- 
particle) Coulomb matrix elements can be expressed as 

V+ = N 4 (4g 2 (l + S 2 )f\ + (1 + g 2 ) 2 F 2 

+ Ag 2 F 3 - 16g 2 F 4 ) , (A3) 
V- =N 4 (l-g 2 ) 2 (F 2 -S 2 F 3 ), (A4) 
U = N 4 ((1 + g 4 + 2g 2 S 2 )F 1 + 2g 2 F 2 

+ 2g 2 S 2 F 3 - 8g 2 F 4 ) , (A5) 
X = N 4 [((1 + g 4 )S 2 + 2g 2 ) F 1 + 2g 2 F 2 

+ 2g 2 S 2 F 3 - 8g 2 F 4 ] , (A6) 
w = N 4 (- 5 (1 + 5 2 )(1 + S 2 )F 1 - 5 (1 + g 2 )F 2 

- 5 (1 + g 2 )S 2 F 3 + (1 + 6g 2 + g 4 )SF 4 ) , (A7) 



with N = \/^J\-2Sg + g 2 and g = (1 - y<l-S 2 )/S. 
Here, we make use of the functions 

Fi = cVb, (A8) 

F 2 = cVbe- bd2 l Q (bd 2 ) , (A9) 

F 3 = cV&e^-^Io (d 2 (b - 1/b)) , (A10) 
F^cVbe'^/^x 

x £ (-l) fe I 2fe ^(26-l/6))l 2fe ^yv / 6 2_ T), 

(All) 

where I„ denotes the Bessel function of n-th order. For 
our purposes, we can neglect terms with \k\ > 1 in the 
sum in F4, since for Tiuiq — 3 meV, B < 30 T, and d = 0.7 
the relative error introduced by doing so is less than 1%. 
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APPENDIX A: HUND-MULLIKEN MATRIX 
ELEMENTS 

Here, we list the explicit expressions for the matrix 
elements defined in Eqs. ^ and ([!(]) as a function of 
the dimensionless inter-dot distance d = a/ae and the 
magnetic compression factor b = + <jJ 2 /u)q where 
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If during the change of v(t) the total spin remains con- 
served, no transitions between the instantaneous singlet 
and triplet eigenstates can be induced during the switch- 
ing. Thus, the singlet and triplet states evolve indepen- 
dently of each other, and the condition on adiabatic 
switching involves Ae (instead of J), i.e. we only need 
to require that l/r s ~ \v/v\ <C Ae/h, which would be 
less restrictive. Also, only f Ts dtj(t) and not J(i) itself 
is needed for the gate operation. Therefore, the adiabatic- 
ity criterion given in the text, while being sufficient, need 
not be really necessary. However, the complete analysis 
of the time-dependent problem in terms of variational 
wave functions is beyond the scope of the present work 
and will be addressed elsewhere. 

We note that it is sufficient to have single-qubit rota- 
tions about any two orthogonal axes. A preferable choice 
here are two orthogonal in-plane axes because magnetic 
fields By parallel to the 2DEG do not affect the exchange 
coupling J(B±) (assuming that we can exclude subband 
mixing induced by a sufficiently strong B\\). 
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